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Abstract 

A study of the particle multiplicity between jets with large rapidity separa- 
tion has been performed using the D0 detector at the Fermilab Tevatron pp 
Collider operating at ^fs = 1.8 TeV. A significant excess of low-multiplicity 
events is observed above the expectation for color-exchange processes. The 
measured fractional excess is 1.07±0.10(stat)loj|(syst)%, which is consistent 
with a strongly-interacting color-singlet (colorless) exchange process and can- 
not be explained by electroweak exchange alone. A lower limit of 0.80% (95% 
C.L.) is obtained on the fraction of dijet events with color-singlet exchange, 
independent of the rapidity gap survival probability. 
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Jet production in hadron-hadron collisions is typically associated with the exchange of a 
quark or gluon between interacting partons. In addition to these color-exchange processes, 
the exchange of an electroweak color singlet (photon, W or Z boson) or a strongly-interacting 
color singlet (such as two gluons in a colorless state), can also produce jets. Two jets sepa- 
rated by a rapidity gap, defined as a region of rapidity containing no final-state particles, has 
been proposed as a signature for jet production via the exchange of a color-singlet (colorless) 
object Theoretical calculations have confirmed that gluon radiation between scattered 
partons is highly suppressed for color-singlet exchange relative to color exchange Hf||. 

The multiplicity of final-state particles in the rapidity interval between jets offers a conve- 
nient way to distinguish color-singlet exchange from color exchange. Color-singlet exchange 
is expected to give a multiplicity near zero for events that have no spectator interactions 
and a minimum bias-like multiplicity distribution for events that contain spectator inter- 
actions In contrast, color-exchange events are expected to have a much higher mean 
multiplicity and to be described by a negative binomial distribution (NBD) or a sum of 
two NBD's (double NBD) |5|,|S|]. Low-multiplicity color-exchange events, which are a back- 
ground to color-singlet exchange, become suppressed as the rapidity interval between the 
jets increases. An excess of low-multiplicity events with respect to the distribution for color 
exchange would indicate the presence of a color-singlet exchange process. 

The magnitude of any observed excess can be used to distinguish between a strongly- 
interacting or purely electroweak color singlet. The exchange of a two-gluon color-singlet, 
which has been proposed as a model for the pomeron, is roughly estimated to account for 
10% of the dijet cross section p],[7|||, while the contribution from electroweak exchange is 
calculated to be about 0.1% 0. The survival probability (S) for rapidity gaps to contain no 
particles from spectator interactions is estimated to be 10-30% @,||,(|. Thus, the fraction 
of dijet events with an observable rapidity gap is expected to be about 1-3% for two-gluon 
color-singlet exchange and 0.01-0.03% for electroweak exchange. 

Although evidence exists for color-singlet (pomeron) exchange in single-diffractive jet 



events RlOlJllI , these events are typically produced with low momentum transfer, 0< \t\ <2 



(GeV/c) 2 . The D0 and CDF Collaborations have previously reported the observation of 



rapidity gaps in dijet events [p~2| , p~3[1 , which, in contrast to single-diffractive jet events, have 
a rapidity gap between the jets and \t\ >900 (GeV/c) 2 . This letter presents a new analysis, 



based on the same data as in Ref. ||12|| , that provides clear experimental evidence for jet 
production via color-singlet exchange at a level inconsistent with an electroweak exchange 
process alone. In this work, pseudorapidity, 77 = — lntan(#/2), is used as an approximation 
for true rapidity. 

The D0 detector and trigger system are described elsewhere [Q. The primary trigger 



used in the present analysis was previously discussed in Ref. |12[ and required "opposite- 
side" jets with a large pseudorapidity separation. In the offline analysis, the two leading E T 
(highest transverse energy) jets are required to have E? > 30 GeV, \r]\ > 2, and f]i ■ f]2 < 0. 
A cone algorithm with radius 1Z = (Ar] 2 + A0 2 )^ = 0.7 is used for jet finding. Events with 
more than one proton-antiproton interaction are removed since extra interactions would ob- 
scure a color-singlet signature and alter the multiplicity distribution. This single interaction 
requirement yields a final data sample of 22,400 opposite-side jet events. 

An independent trigger required events on the "same-side" of the detector in pseudo- 
rapidity. A final sample of 23,200 same-side jet events is obtained by requiring a single 
interaction and two jets with Et > 30 GeV, |?7| > 2, and 771 • 772 > 0. The same-side sample 
provides a qualitative measure of the color-exchange background multiplicity in the central 
rapidity region due to the color flow between the scattered and spectator partons. Hard sin- 
gle diffraction, which could produce a central rapidity gap with two forward jets, is highly 
suppressed by the trigger which required a coincidence of hits between the forward and 
backward luminosity counters (1.9 < <4.3). 

The electromagnetic (EM) section of the calorimeter is used as the primary means of 
measuring the particle multiplicity. The EM calorimeter has a low level of noise and the 
ability to detect (with an energy-dependent efficiency) both neutral and charged particles 
for 1 77 1 < 4. A particle is tagged by the deposition of more than 200 MeV transverse energy in 
an EM calorimeter tower (Ai]xA<j) = 0.1x0.1). The central drift chamber (CDC) is efficient 



for detecting charged particles for <1.3 and provides an independent measurement of the 
multiplicity. 

Figure [I] shows the number of EM calorimeter towers above threshold (n ca i) versus the 
number of CDC tracks (n^k) i n the region \rj\ < 1.3 for the (a) opposite-side and (b) same- 
side samples. The two distributions are similar in shape except at very low multiplicities, 
where the opposite-side sample has a striking excess of events, consistent with a color-singlet 
exchange process. For both samples, n ca i and n tr k are strongly correlated, confirming that 
either can be used as a measure of particle multiplicity. 

A model for the multiplicity in color-exchange events is necessary to measure the low- 
multiplicity excess observed in the opposite-side data. We use the double NBD, which has 
four parameters and a relative normalization, to parametrize the color-exchange component 
of the opposite-side multiplicity distribution between jets. 

This parametrization is supported by color-exchange data and Monte Carlo samples. 
Figure 0(a) shows that the same-side n ca \ distribution for \rj\ < 1 is well-parametrized by 
the double NBD over the full range of multiplicity giving a x 2 /df (degree of freedom) = 0.9. 
Additional support is given by a color-singlet-exchange-suppressed subset of the opposite- 
side data obtained by demanding the presence of a third jet (with E T >8 GeV) between the 
two leading jets. Figure ||(b) shows the multiplicity for these events in the pseudorapidity 
region between the jet cone edges, excluding the multiplicity of the third jet. The double 
NBD fits the distribution reasonably well {x 2 /df = 1.2) over the full range of multiplicity. 
Monte Carlo color-exchange events, generated with HERWIG [[15] and PYTHIA |Tj| and passed 



through a simulation of the D0 detector |5|,[17|, provide further support for the double NBD 
(not shown). None of these color-exchange samples have a significant excess of events at low 
multiplicity from physics processes or detector effects. 

The low- multiplicity excess observed in Fig. |I|(a) is determined as follows. First, a double 
NBD is fit to the opposite-side calorimeter multiplicity distribution between the two leading 
jets using a binned maximum likelihood method. The overall normalization is set to the 
number of events in the fitted region and the relative normalization of the two NBD's is 
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determined by minimizing a modified Kolmogorov-Smirnov (K-S) statistic, which is defined 
as the sum of the squared differences between the cumulative distributions of the data and 
the fit. Next, the starting bin (no) for the fit is incremented successively by one and the 
remaining bins refit until the x 2 jdf i n the first five bins of the fit is less than 2.0 (i.e. 
C.L. > 7%). The final fit is then extrapolated to zero multiplicity to determine the excess 
above the expected color-exchange background for n ca i < no. Note that any experimental 
effects that produce a smearing of zero multiplicity color-singlet events to slightly higher 
multiplicities are reduced by the integration over low-multiplicity bins. 

Figure |3](a) shows the n c& \ distribution between the cone edges of the two leading jets for 
the opposite-side data sample. Also shown is the double NBD fit (x 2 /df = 0.9, n = 3) and 
its extrapolation to zero multiplicity. In contrast with the color-exchange data and Monte 
Carlo samples, the opposite-side sample has a significant excess at low multiplicity. The mul- 
tiplicity in the central calorimeter (\n\ < 1), a region well away from the jet edges, is plotted 
in Fig. 0(b) for the same events. Although the distribution has a lower mean multiplicity, 
its shape (excluding low-multiplicity bins) is also well-fit by a double NBD (x 2 /df = 1.0, 
no = 2), and a clear excess is again present at low multiplicity. The excess, defined as the 
integrated difference between the data and the double NBD fit in the extrapolated region, 
is 225 ± 20 and 237 ± 21 events for Figs. ||](a) and (b), respectively. 

Using a single NBD instead of a double NBD to fit the data in Figs. |3|(a) and (b), as done 



in Ref. JT8[, gives a significantly larger excess of 339 and 421 events, respectively. Although 
the single NBD fit gives a x 2 /d/=l.l, K-S tests show that the shape of the double NBD is 
significantly favored over the single NBD. 

The excess observed at low multiplicity is consistent with the presence of color-singlet 
exchange, but the fractional excess is a more relevant quantity for theoretical comparisons. 
The fractional excess f$ is defined as the average excess above the double NBD fit for starting 
points of no and uq+1 divided by the total number of events. We obtain a fractional excess 
of fs = 1.07% from the multiplicity distribution between the jet edges shown in Fig. |3](a), 
and a value of 1.14% using the distribution in the region |?7| <1 from Fig. ||](b). 
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The exact fit parameters depend on the definition of particle multiplicity, but the frac- 
tional excess is relatively independent of this definition. Varying the calorimeter tower 
threshold and refitting the multiplicity yields a consistent value of f$. Redefining a particle 
as a CDC track or a "cluster" of neighboring calorimeter towers results in an excess that 
is 3-15% greater. We thus assign a conservative +15% systematic error on f s due to the 
uncertainty in defining particle multiplicity. 

Systematic effects from uncertainties in jet reconstruction, energy scale, and acceptance 
result in a ±5% uncertainty in the measured excess. Uncertainty in the identification of 
single-interaction events gives an additional error of 1l\ 4 %. A systematic error of ±10% 
due to the uncertainty in the background is determined by varying the fit parameters to 
produce a one unit change in y 2 an d remeasuring the excess. Combining the above errors 
in quadrature gives / 5 = 1.07±0.10(stat)ig;!(syst)%. 

Although a low-multiplicity excess is clearly present in the opposite-side sample, theoreti- 
cal interpretation of fs is complicated by uncertainties in the survival probability. In contrast 
with our previous measurement which used the number of events with zero multiplicity to 
place an upper limit on the rapidity gap fraction ((cr singlet /(x) xS< 1.1% at 95% C.L.) JT2| , 
the measurement of fs may include some portion of the color-singlet events with low- 
multiplicity spectator interactions. Nevertheless, the value of fs is consistent with the 1-3% 
expected for strongly-interacting color-singlet exchange. The measured value of fs is used to 
obtain a lower limit of 0.80% (95% C.L.) on cr sing i e t/a, the fraction of dijet events produced 
via color-singlet exchange, independent of the actual value of the survival probability. 

The opposite-side data can also be used to exclude pure electroweak exchange combined 
with color-exchange background as the source of the excess. A PYTHIA Monte Carlo study 
using simulated D0 jet acceptance and efficiency gives a value of 0.09% for the fraction of 
dijet events with electroweak exchange, which is comparable to the result in Ref. |7| but 
includes higher-order radiative corrections. We assume a survival probability of S = 100% 
to determine the maximum expected number of events from electroweak exchange, and add 
this to the number of color-exchange background events from the fit. The statistical and 
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systematic uncertainty in the expected number of events is then used to determine the prob- 
ability for a fluctuation to the number of observed events. Using only the zero multiplicity 
bin, the probability that the observed excess in Fig. |3](a) is due to the combination of elec- 
troweak and color exchange is less than 10~ 10 . Adding other low-multiplicity bins further 
decreases this probability. 

In conclusion, we have presented evidence for strongly-interacting color-singlet exchange 
from a study of opposite-side dijet events with jet Et > 30 GeV and \r]\ > 2. A strik- 
ing enhancement of low-multiplicity events is observed independent of the details of the 
color-exchange background parametrization. The double NBD parametrization provides a 
measurement of the color-exchange background and confirms that this background is small 
for the rapidity gap fraction previously measured in Ref. ||12|| . The fractional excess above 
the color-exchange background is found to be 1.07 ± 0.10(stat)l ;i3( s y s t)% which is consis- 
tent with the presence of strongly-interacting color-singlet exchange and inconsistent with 
electro weak exchange alone. 
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FIGURES 




FIG. 1. The calorimeter tower multiplicity (n ca i) versus the charged track multiplicity (n tr k) in 
the pseudorapidity region \rj\ < 1.3 for the (a) opposite-side and (b) same-side samples as described 
in the text. 
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FIG. 2. The calorimeter tower multiplicity (n ca i) for the (a) same-side and (b) color-singlet 

exchange suppressed data samples. The solid lines are double NBD fits to the data. The insets 

show each plot on a log-log scale. 
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FIG. 3. The number of events versus n ca i (shifted up by half a unit in multiplicity) is plotted 

on a log-log scale for all opposite-side jet events, where (a) shows the multiplicity between the 

cone edges of the two leading jets and (b) shows the multiplicity in the region \rj\ < 1. The solid 

lines represent the double NBD fits for (a) no = 3 and (b) uq = 2, while the dashed lines show the 

extrapolation of each fit to n ca \ = 0. 
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